Lieb-Schultz-Mattis theorem for quasi-topological systems 
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In this paper we address the question of the existence of a spectral gap in a class of local Hamiltonians. These 
Hamiltonians have the following properties: their ground states are known exactly; all equal-time correlation 
functions of local operators are short-ranged; and correlation functions of certain non-local operators are critical. 
A variational argument shows gaplessness with CO <>= fc^ at critical points defined by the absence of certain terms 
in the Hamiltonian, which is remarkable because equal-time correlation functions of local operators remain 
short-ranged. We call such critical points, in which spatial and temporal scaling are radically different, quasi- 
topological. When these terms are present in the Hamiltonian, the models are in gapped topological phases 
which are of special interest in the context of topological quantum computation. 
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I. INTRODUCTION 

Recently, there has been great interest in various uncon- 
ventional states of matter, in particular those that might ap- 
pear in strongly-interacting 2D systems. These exotic phases 
may be characterized by unconventional order parameters. 
Alternatively, the phases may be topological. In this case 
there is no local order parameter, but there may be a non- 
local order parameter, which is related to the topological 
properties of the state: exotic braiding statistics, quantum 
number fractionalization, and a ground- state (GS) degener- 
acy which depends only on the topology of the underlying 
2D manifold^. The low-energy effective description of these 
phases is a topological quantum field theory (TQFT). A vari- 
ety of topological phases are known to exist in the quantum 
Hall systems^. It has been conjectured that such phases also 
occur in frustrated magnets, where they may be connected 
to superconductivit y^''*'^'^'^-^-^ . Topological phases are also 
an attractive platform for quantum computation, where their 
insensitivity to local perturbations leads to fault-tolerance'>S. 
Non-Abelian topological phases are particularly interesting in 
this context because, in many of these phases, the braiding 
of quasiparticles generates a set of transformations which is 
sufficient for universal quantum computationii. In this pa- 
per, we shall concentrate on microscopic models which are 
related to a class of P and invariant non-Abelian topo- 
logical phasesi^. 

A set of conditions which places a microscopic model in 
such a topological phase can be briefly summarized as fol- 
lows. We suppose that the low-energy Hilbert space can be 
mapped onto that of a quantum loop gas. This is the case 
in a large class of models, including dimer models, certain 
spin models, and some interacting hard-core boson models. 
In such a model, basis states are associated to collections 
of non-intersecting loop a'^''^''^ . We give some examples in 
section |III] A Hamiltonian can act on states in this Hilbert 
space by doing the following: (i) the loops can be contin- 
uously deformed - we will call this an isotopy move; (ii) 
a small loop can be created or annihilated - the combined 



effect of this move and the isotopy move has been dubbed 
"c/-isotopy"i2ii^ii^; (iii) finally, when exactly k+ \ strands 
come together in some local neighborhood, the Hamiltonian 
can cut them and reconnect the resulting "loose ends" pair- 
wise so that the newly-formed loops are still non-intersecting. 
More specifically, in order for this model to be in a topological 
phase, the ground state of this Hamiltonian should be a super- 
position of all such pictures with the additional requirements 
that (i) if two pictures can be continuously deformed into each 
other, they enter the GS superposition with the same weight; 
(ii) the amplitude of a picture with an additional loop is d 
times that of a picture without such loop; (iii) this superposi- 
tion is annihilated by the application of the Jones- Wenzl (JW) 
projector that acts locally by reconnecting k+\ strands (for a 
detailed description see the following section). The main goal 
of this paper is to investigate the energy spectrum of a system 
subject to the first two conditions. 

We shall show that a generic local Hamiltonian which en- 
forces t/-isotopy for its ground state(s) is necessarily gapless 
provided that \d\ < \/2. (We shall also argue that, surpris- 
ingly, even the addition of the JW projector to the Hamilto- 
nian will not open a gap for \d\ — %/2.) Such a Hamiltonian 
is a sum of projection operators (enforcing on every plaque- 
tte of the lattice both isotopy invariance and the value d for 
a contractible loop). These projection operators do not com- 
mute with each other, but they are compatible with each other 
in the sense that they all annihilate the ground state. (Such 
Hamiltonians have also arisen in the context of the quantum 
Hall effect, where the Haldane pseudopotentials are projec- 
tion operators which annihilate the Laughlin states'^-'\ but do 
not commute so the excited states are not known exactly, and, 
under a general name of "parent Hamiltonians", in quantum 
antiferromagnetisroi^ii^i^.) Exact knowledge of the ground 
state enables us to construct a variational ansatz for the low- 
est energy excited state. The strategy is quite similar to the 
single-mode approximation (SMA), in which pf/|0) is the trial 
excited state, where p is some conserved charge. In the case 
of a broken symmetry state, p is the charge which generates 
the symmetry transformation which is spontaneously broken 
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in the ground state. Our method generalizes the SMA in an 
important way: we do not rely on the existence of any con- 
served charges and, indeed, the model need not have any. In- 
stead, we have the less restrictive condition that the configura- 
tion space of the model should break into two (or more) sec- 
tors whose volume is parametrically larger than the boundary 
between them, e.g. by a factor of the system size L. If the 
Hamiltonian only has matrix elements between nearby points 
in configuration space, then (a la Lieb-Schultz-Mattis^-) we 
can construct a wavefunction which is equal to the ground 
state wavefunction except for a relative sign change between 
the two sectors. The energy cost of this "twisted" excited state 
would be at most ~ 1 /L and is even smaller in the case of the 
models which we consider here. On the other hand, if the 
Hamiltonian directly connects the two putative sectors, then 
it means that in reaUty these are not two distinct sectors, and 
a sign change in the wavefunction would have an energetic 
penalty which does not scale to zero as L — * oo. 

In this paper, we prove the general result on gaplessness 
described in the preceding paragraph. We apply this result 
to microscopic models implementing li-isotopy and find the 
remarkable result that they are gapless with (i>°^ in spite of 
the absence of any power-law equal-time correlation functions 
of local operators. 



n. rf-ISOTOPY AND ITS LOCAL SUBSPACES 

In this section we present a more formal definition of d- 
isotopy ; it can be skipped on a first rea ding. R eaders interested 
in the details are also referred to Refs. Il2il3lfl5[ 

Although we will eventually be considering a system on a 
lattice, it is useful to begin by defining the Hilbert spaces of 
interest, and Vd, in the smooth, lattice-free setting. Con- 
sider a compact surface Y and the set S of all multiloops^ 
X <zY . If dY (the boundary of Y) is non-empty, we fix once 
and for all a finite set P of points on dY with X{^dY =P. We 
assume Y is oriented but X should not be. There is a large vec- 
tor space C'^, of complex valued functions on S. We say X and 
X' are isotopic (X^X') if one may be gradually deformed into 
the other with, of course, the deformation being the identity 
on dY (see Fig.[T]). 




FIG. 1: Isotopy on an annulus: X ^ X' but X X" . 

We may view the isotopy relation as a family of linear con- 
straints on C'^, namely =0 if X-X'. The sub- 
space satisfying these linear constraints is now only of count- 
able dimension; it consists of those functions which depend 
only of the isotopy class [X] of X and may be identified with 
where [S] is the isotopy classes of multiloops with fixed 
boundary conditions. Note that all isotopes can be made by 



composing small locally-supported ones, so the relations we 
just imposed are "local" in the sense that they can be imple- 
mented with purely local terms in the Hamiltonian. 

Let us go further and define an additional local relation 
which, when added to isotopy, constitutes the "t/-isotopy" re- 
lation. This relation is: 



c/'P(X)-^(XUO) =0 



(1) 



It says that if two multiloops are identical except for the 
presence of a small (or, it follows, any contractible) circle, 
then their function values differ by a factor of d, a fixed pos- 
itive real number. In cases of interest to us \ < d < 2, s,o 
our function is either neutral to or "likes" small circles. We 
call the subspace obeying all these constraints the ^f-isotopy 
space of Y (with fixed boundary conditions) and write it as 
Vd C Ct'^l C C-^. The subspace Fj(r2) is still of countable 
dimension, or extensively degenerate, on the torus T~. The 
vector space Vd is clearly the ground state manifold (GSM) 
of a local Hamiltonian acting on Cl"^! . 

It is a remarkable facti^ill that it is very difficult to add 
any further local relations to J— isotopy without killing the 
vector space entirely. For the physically interesting cases a = 
g7ci/(<r+2)^ A; = 1,2,3,..., there is such a local relation and a 
natural positive definite inner product on V d (see refs. [l2lfl5l 
for definition). 

In these cases the local relations are essentially the Jones- 
Wenzl idempotents: 
k=\: 
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0, (4) 



see Ref.|22|for a recursive formula. These relations define a 
finite dimensional Hilbert space Vd{Y) C Vd{Y) C Cl'^l C C'^. 

In Refs. lTWlsl it is explained that VdiY) is the Hilbert spaces 
for doubled SU{2)k Chern-Simons theory on Y. It has been 



3 



arguedii that a Hamiltonian with a ground state manifold cor- 
responding to Vci is perched at a phase transition. When per- 
turbed, (infinitesimally for A: = 1 or 2 and under a larger defor- 
mation for A: > 3) it will go into a topological phase with low- 
energy Hilbert space Vd, i.e. into a phase described by doubled 
SU{2)k Chern-Simons theory. Very briefly, the Hilbert space 
of a TQFT such as doubled SU{2)k Chern-Simons theory, can 
always be defined as the joint null space of commuting lo- 
cal projectors^, implying the existence of a local Hamilto- 
nian with a spectral gap in the thermodynamic limit. Once a 
Hamiltonian H^i has imposed (i— isotopy, i.e. GSM{Hj) = Vd, 
an extensive degeneracy has been created; the only local way 
of lifting this extensive degeneracy (to a finite degeneracy) 
without creating frustration^ is to add the Jones-Wenzl pro- 
jector to such a Hamiltonian. 



III. MICROSCOPIC LATTICE MODELS 

We now briefly review some examples of microscopic 
Hamiltonians whose ground state(s) are described by d- 
isotopy. These Hamiltonians may not be particularly re- 
alistic (e.g. the two spin Hamiltonians presented here do 
not conserve the total spin), but they have the advantage of 
being relatively simple local Hamiltonians with the desired 
ground states. These ground states have a nice property that 
their square-norms are in one-to-one correspondence to the 
partition functions of known statistical mechanical models 
- a so called plasma analog-^^ This approach has proved 
useful in studying many frustrated quantum models and, 
in particular, those with topological and quasi-topological 

The first example was presented in Ref.[l3 and was inspired 
by Kitaev's modeU^. The model is defined on a honeycomb 
lattice, with the elementary degrees of freedom being s ~\/2 
spins situated on its links (alternatively, one can think of these 
spins as occupying the sites of a kagome lattice, but the former 
description lends itself nicer to a loop representation). The 
Hamiltonian is given by 



H 



1) 



o^t 



d P 



-Fl{Fl)'-Fl-{Fp 



d 

2\t 



-fUfI) 



(5) 



where 9{, (v) is the set of 3 links neighboring vertex v, and 



Fp = OiO^o^a^a.a^ 

Fp ~ Oj^Oj 0304 a^Og + cyclic perm. 



- cyclic perm, 
cyclic perm. 



Here, 1 , 2, . . . , 6 label the six edges of plaquette p. This Hamil- 
tonian is a sum of projection operators with positive coeffi- 
cients and, therefore, is positive-definite. The eigenvalues of 
the first term in (|5]l are 2,0, corresponding, respectively, to 
whether there is an even or odd number of o~ = — 1 spins 
neighboring this vertex. When eigenvalue zero is obtained 
at every vertex, the = 1 links form loops. Hence, the zero- 
energy subspace of the first term is spanned by all configu- 
rations of multiloops (on the honeycomb lattice, they cannot 
cross). 

The term on the second line of (|5]) is a projection operator 
which annihilates a state |^) if the amplitude for all of the 
spins on a given plaquette to be up is a factor of d times the 
amplitude for them to all be down, i.e. if the amplitude for a 
configuration with a small loop encircling a single plaquette 
is a factor of d times the amplitude for an otherwise identical 
configuration without the small loop. The other three Unes of 
the Hamiltonian vanish on a state \^) if it accords the same 
value to a configuration if a loop is deformed to enclose an 
additional plaquette - in other words, these terms enforce the 
usual isotopy relations. These terms are graphically repre- 
sented in Fig.|2l 



: d 




(6) 



FIG. 2: The action of the terms in the last four lines of the Hamilto- 
nian l|5j represented graphically. The solid bonds correspond to the 
up-spins. Notice that the application of the Hamiltonian to any of the 
above plaquette configurations results in a superposition of this con- 
figuration and its counterpart with the appropriate amplitudes. The 
relative amplitudes of these configurations in the ground (0 eigen- 
value) state correspond to the "d"-isotopy rules described in the In- 
troduction. 

The second example was presented in Ref. [isi motivated 
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by the connection to a classical statistical mechanical model, 
namely the self-dual Potts model (this connection will be ex- 
plored in the next section). The model is defined on a square 
lattice; the degrees of freedom are once again i = 1/2 spins 
situated on its links. The Hamiltonian is given by 




in the notation of Ref. [HI The sum in Eq. ^ is taken over all 
elementary plaquettes of the lattice, |3) is a state with (any) 
three up-spins around a given plaquette, |4) has all four spins 
up. The second sum is taken over all vertices; |1) and |0) 
correspond to a single spin or no spins up around a given ver- 
tex. Notice that this Hamiltonian is self-dual under flipping 
all spins and going to the dual lattice. 

The loops are now defined on the surrounding (or midpoint) 
lattice, i.e. the lattice obtained by connecting the mid-points 
of adjacent edges. One should think of placing a double-sided 
mirror along a bond whose spin is up and placing a mirror 
along a dual bond if the spin is down. Loops are formed by 
propagating light in this labyrinth of mirrors. The action of 
the first term Eq. (|7| is demonstrated in Fig. [3] The action of 
the second (vertex) term is completely analogous due to the 
aforementioned duality. 




FIG. 3: The action of the first (plaquette) term of Hamiltonian l|7j on 
a shaded plaquette. A small loop is created or annihilated (with the 
appropriate amplitude) by being merged with another loop. 

Finally, a reader interested in a more realistic Hamiltonian 
is referred to |1T3,3 1] where a possibility of finding a J-isotopy 
ground state(s) in an extended Hubbard model is discussed. 
We shall refrain from reviewing that construction here due to 
its complexity, unnecessary for the purpose of this paper 

Let us turn to the common features of the above Hamilto- 
nians. Both the terms on the final four lines of (|5]l and the 
plaquette and vertex terms of ^ do not commute with each 
other However, they are compatible in the sense that each of 
these terms annihilates the following state 

i>i'o)-i:^''^"V)- (8) 

{«} 

which is a superposition of all configurations a of multiloops 
weighted by a factor of d to the number of loops £{a). 

Notice one important difference between these examples; 
by construction, in the second case the loops are fully packed: 
every edge of a surrounding lattice is traversed by a loop. Be- 
cause of this constraint, the notion of simple isotopy is mean- 



ingless here; the multiloops, however, can fluctuate by "eject- 
ing" and "absorbing" small loops. On an L x L torus, the 
ground state degeneracy is ^ because the Hamiltonian does 
not mix states |a) with different winding numbers^^. The dif- 
ferent ground states are given by ^ but with the sum over a 
restricted to a single topological class. Notice that at least in 
the first example, the Hamiltonian dynamics is ergodic within 
a given topological sector; this is particularly trivial to see 
for the zero winding case since every loop configuration can 
be reduced to a configuration with no loops by applying the 
moves depicted in Fig. |2]to first shrink the loops to a single 
plaquette and then to annihilate them. 



In the second model, ergodicity is less obvious since the 
loop model is fully-packed and thus reducing any given con- 
figuration to a state with no loops is impossible. Moreover, it 
has been observed in Ref. [Tsl that for a non-zero winding num- 
ber for certain finite tori there are configurations that are not 
connected by the the "moves" of the Hamiltonian ([7]). For the 
zero winding number, however, it can be shown that all config- 
urations are connected to the state with all spins down, which 
in turn translates into a maximum possible number of loops, 
one per each plaquette of a dual lattice. Hence the Hamilto- 
nian is ergodic in the zero winding sector, and this is the sector 
we will be concerned with for the remainder of the paper (Er- 
godicity is important because the twisted states we propose 
should have nothing to do with degeneracies associated with 
possible non-ergodicity.) 



More generally, the ground state on any genus g > I sur- 
face is infinitely degenerate in the thermodynamic limit. As 
we saw in the previous section if d ^ ±1 oi d = ±V2, there 
is a Jones-Wenzl projector which also annihilates the ground 
state ([8]) on a topologically-trivial manifold but mixes differ- 
ent winding number sectors on higher-genus surfaces. Hence, 
in either model, at d ~ ±l,±\/2, there are two Hamiltonians 
which have the same ground state on the sphere (or, equiva- 
lently, in the zero winding sector). The first ones, given by 
either Eq. ^ or Eq. (|7|i, have extensively degenerate ground 
states on the torus while the other kind (with the Jones-Wenzl 
projector added) have finite degeneracjii^ The second type 
leads to a topological phase with an energy gap for d — 1—, 
where the resulting model is exactly soluble, and, we believe, 
for d = —1 (the situation at li = ±\/2 is more subtle and will 
be addressed later). The spectrum of the first kind is the main 
subject of this paper. 



While the ground state can be obtained exactly, excited 
states cannot because the different operators in both (|5]l and 
(|7]i do not commute with each other We will use a varia- 
tional ansatz to show that in the absence of other terms, such 
as Jones-Wenzl projectors, the spectra of these Hamiltonians 
are gapless in the zero winding sector. But before we can pro- 
ceed, we need to address the inherent structure of these ground 
states in more details. 
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IV. MAPPING OF THE GROUND-STATE TO A 
STATISTICAL MECHANICS PROBLEM 

Many properties of the ground-state wavefunction can be 
obtained by observing that the norm of the ground state is 
equal to the partition function of a classical loop model: 

{a} 

where a denotes a particular configuration of loops on a lattice 
(a "snapshot" of multiloops). The specific details of possible 
configurations depend on a particular choice of a Hamiltonian; 
in what follows we shall consider the cases relevant to each of 
the proposed Hamiltonians. 



A. Potts model: random clusters and loops 

The ^-state Potts model, originally introduced as a general- 
ization of the Ising model, is defined by the following classical 
Hamiltonian: 

-p^=/^8a„a,, (10) 
< '■,;■> 

where the sum is taken over all pairs of nearest neighbors and 
O/ is a discrete "spin" variable that can take on q different 
values, e.g. O/ = 1 , . . . , <7. 

We now review the the Fortuin-Kasteleyn {AKA random 
cluster) representation for the ^-state Potts model and under- 
line its basic properties. 

Given the Hamiltonian [TO] the partition function can then 
be written as 

Zpo„. = Ee-P-'^=^ n e''"' "^ 

{o} {o} <'■../■> 

[l + (e^-l)5a,.a,] 

{o} <i j> 

=L n [^+^'K.o^ : (11) 

{o} <ij> 

with V = e'' — 1 . 

The next step is to expand the product in Eq. ( fTTI) . Every 
pair of nearest neighbors contributes a factor of either 1 or v to 
each of the resulting terms, with the latter possibility available 
only if the neighboring spins agree (a, = Oj). Therefore every 
such term has a simple graphical representation: a bond is oc- 
cupied if the pair of spins it connects contributes a factor of v 
into a given term, and it is left vacant otherwise. Since a bond 
can be placed between the sites only if their spins "agree", all 
spins belonging to the same bond cluster must have the same 
value. 

The partition function thus becomes 

Zpotts = E n [1 +^'5a,.a,] = E L v*(")A(a,(o) , (12) 

{a} <i,j> {o} {to} 



where b{({)) is the total number of occupied bonds in a bond 
configuration co, and A(a,a)) is the appropriate collection of 
Kronecker delta-symbols that enforces the "agreement" be- 
tween the spin and the bond configurations. 

The next step is to change the order of summation in 
Eq. ( fT2b (which is fine as long as the lattice is finite) and then 
sum over all spin configurations a. The only constraint on the 
spin variables (for a fixed bond configuration co) is the one 
that has been mentioned earlier: all connected spins must take 
on the same value (one of the q possible). Therefore every 
connected cluster, as well as each isolated site, contributes a 
factor of q to the resulting bond weight, and the partition func- 
tion becomes 

Zpotts = Ev^(»)^^("), (13) 
{to} 

where c(a)) is the total number of connected components (in- 
cluding isolated sites). The partition function ( fTSl l defines FK 
(random cluster) representation^^^^ 

While all previous considerations applied to any number 
of spatial dimensions, in 2D the partition function (fTSI l can 
also be written in a way which makes its self-dual property 
apparent. In order to do this, notice that if we think of an 
unoccupied bond of the actual lattice as an occupied bond of 
the dual lattice (therefore b* =B — b), then every circuit (face) 
of the actual lattice contains a dual connected component: / = 
c* (here we define b* = b{&*), c* = c(co*) with CO* denoting 
the configuration of dual bonds). Therefore, using the Euler 
relation (for a planar graph) 

=b{(£,)+c{&)-N (14) 

with /(co) being the number of circuits (defined as a minimum 
number of bonds that one has to cut in order to make a graph 
consist only of tree^) and being the total number of sites, 
we can rewrite Eq. ( [13] ) as 

zpotts = I vf-^^'q^ = E v-f {^y = v^' {^y 

{tn} {to} ^' {to} ^' 

(15) 

If we ignore the uninteresting analytic prefactor in Eq. ( fTSl ). 
we immediately see that the system is self-dual when v ^q/v, 
i.e. V = 

A so-called polygon decomposition^ lets us relate random 
clusters to a loop gas on the surrounding lattice. We think 
of an occupied bond as a double-sided mirror placed at the 
site of the surrounding lattice. If a bond is not occupied, then 
its dual bond is considered a mirror. Thus every site of the 
surrounding lattice gets one of the two possible mirrors. We 
then use these mirrors to construct paths as shown in figured 
Since these paths have no sources or sinks, they always form 
loops that either surround the clusters or are contained inside 
clusters (in the latter case, the loops surround dual clusters). 
The number of loops I is then given by £ = c + c* . If v = ^ 
- i.e. if the Potts model is at its self-dual point, then 

Zself-Dual = E(V^)'^"*- d^) 

{«} 
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FIG. 4: A typical cluster configuration for the Potts model is shown 
by dashed lines. Spins belonging to the same cluster take the same 
value, which must be summed over the q possible values, as de- 
scribed in the text. Clusters can be represented by loops on the sur- 
rounding lattice, shown by solid lines. 

where the sum is taken over all fully-packed loop configura- 
tions on the surrounding lattice. Clearly, with the choice of 
q = li"*, the partition function (fT6] l becomes the norm of the 

(2) 

ground State (^o|^o) of the quantum Hamiltonian //^ . 

B. Random clusters on a torus 

In this section we will take the reader through some mathe- 
matical details in order establish certain properties of the crit- 
ical FK model. These properties will be relied upon later, in 
the course of demonstrating the central result of this paper 

For the sake of concreteness, let us focus the FK representa- 
tion of the critical (7— state Potts model (1 < < 4) on a square 
lattice with periodic boundary conditions in both directions, 
i.e. the torus. In the case of q~ 1, the statistical mechanical 
model reduces to critical bond percolation. In this section we 
are focusing on the FK clusters established in Sec. lIV Al 

Let us begin by proving: 

Proposition 1. Fixing 1 < ^ < 4 and A, > there is an £> 
so that for all L sufficiently large there is a probability greater 
than e that the largest FKq cluster in the Lx L torus has Eu- 
clidian diameter < XL. 

Proof. The technology for this type of result was discov- 
ered and developed by the 1990's in the context of criti- 
cal percolatio n^^-^^-^^'^^ but it can be extended to other crit- 
ical systems provided analogs of the Russo-Seymour-Welsh 
(RSW) inequality on crossing probabilities for rectangles and 
of the Fortuin-Kasteleyn-Ginibre (FKG) inequality ("mono- 
tone events are positively correlated") hold. The proof strat- 
egy is to build some wiggly approximation Gw, lying entirely 
in the complement of FK clusters, to a rectilinear grid of scale 
(roughly) Gw must not be too wiggly as we need the boxes 
in its complement to have diameter < X (see Fig. [S]). Since 
each FK cluster lies in some such box, the clusters also have 
diameter < X. 

Building the grid Gw is a random construction and must 
succeed with probability > e. The naive idea for building the 




FIG. 5: Construction of a wiggly approximation Gw to a rectangular 
grid as described in the text. Gw is built as a subset of the dual clus- 
ters; the goal of the construction is to "trap" the actual FK clusters 
inside Gw, as shown here. 



grid is to first condition on no clusters meeting the (roughly) 

(si)" ^^'•^^ ^^-^^ 4-coordinated vertices) of the desired grid Gg. 
Then, use RSW to further condition on cluster-disjoint arcs 

lying in small rectangles near the (roughly) 2 (^) bonds of 
Gs- What goes wrong is that the individual sites are simply too 
small; in the analogous electromagnetic problem the capaci- 
tance vanishes in the thermodynamic limit. The resolution is 
to consider small annuli enclosing each site and substitute for 
the site condition the event that an essential circle percolates 
around the annulus (also an RSW result). 

Now, elementary planar topology allows bits and pieces of 
the annular and rectangular percolation paths to be hooked to- 
gether to build the desired Gw as shown in Fig.|5] Throughout 
this construction, it is essential that these percolation events 
are at least non-negatively correlated (we remind the reader 
that the properties of dual clusters are identical to those of the 
actual FK clusters at criticality). This is precisely what the 
FKG inequality provide o^^'^^ . □ 
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An immediate consequence of proposition [T| is that there 
is a non-zero probability that an arbitrary multi-loop is in the 
zero-winding-number sector. Thus, we may confine our dis- 
cussion to that sector and use the measure on that sector in- 
duced by the measure c/^^f") on all multi-loops. 

Let us pause for some definitions. Given a loop a in the 
torus r, let I |a| |x, or simply | |a| | for short, be the "breadth in 
the x-direction of a" where a represents a (any) lift of a up to 
the universal cover ^ T^. Recall that the universal cover 
of the torus is constructed by unwrapping completely in 
both the X and y directions. Very concretely, we may build the 
cover from our original LxL square by taking a copy for 
each Gaussian integer and gluing these together (correspond- 
ing to Gaussian integers differing by 1 or /) to form a tiling 
of the plane. (Technically, it is the discrete plane I? that we 
produce.) A loop a in may be regarded as a map from 
the circle 5' into and a lift a is any map to M? making the 
following diagram commute: 



then it descends one-to-one into the torus, implying area(X) < 
1. Thus, 2r8 < 1 and so r < □ 

Now apply the lemma to arcs y within an FK cluster K (in 
the trivial sector) on T . Since we are in the trivial sector, we 
may define \ \K\ \ exactly as we defined the norm on loops. By 
the lemma, if \\K\ \ =8, K will come within distance 1 /8 of 
completing a nontrivial wrapping of the torus. An application 
of the Russo-Seymour- Welsh (RSW) inequality shows that it 
is unlikely (vanishing algebraically in 1 /8) to come close to 
wrapping but fail to wrap completely. RSW amounts to 
integrating the effect of bringing bonds, which may join large 
clusters, in and out of our FK snapshot. Such fluctuations 
cannot be created by any local for since the weight 
of the update is nonlocal. Nevertheless, these fluctuations do 
preserve the measure c/^^'"). So, the fraction of critical, topo- 
logically trivial, FK^ snapshots (1 < q <A) which contain a 
cluster/: with ||/:|| > 8 is ^ijgy. □ 



R2 

~ / 
a / 
/ 

/ 

/ Y 

1 T 

5' > 

Concretely, the image of a is obtained by taking an arc in 
the Lx L square and following its continuation in the tiling 
of until it closes back on itself. The assumption that the 
winding sector is trivial means any such a is finite and in fact 
contains exactly as many bonds as a contained in T^, i.e. < 
2L2. 

To return to the definition of | |fl| |, we say | |fl| | = max — 
Xj\, where the maximum is taken over all sites on a lift 

a. Note that any two lifts a and a are congruent by a translation 
of M^, so ||fl|| is well defined. 

Finally, for a multi-loop a C we define ||a|| = 
max„ca \ We are now ready to state a proposition about 
the diameter of these unwrapped components of random 
multi-loops from the trivial sector. 

Proposition 2. In the scaling limit, ||a|| is almost surely a 
bounded function on the trivial winding sector More pre- 
cisely, the probability that ||a|| > rL has an upper bound 
which is independent ofL. 

Proof. We need a geometric lemma. 

Lemma 1. Let T be a Euclidian torus of area 1. Let y : 
[0, 1] T be an arc so that the (Euclidian) distance between 
lifted endpoints in the universal cover, d(^{0) ^y{\)) = 8. Then 
there is a nontrivial deck translation ( additive action of a 
Gaussian integer) "/ ofyso that dist(y,y') < ^. 

Proof Let X be the radius r neighborhood of y in the in the 
cover M?. Integrating the lengths of slices of X perpendicular 
to the straight line segment [y(0),y(l)] we see (by Fubini's 
theorem) that area(X) > 2r8 (where area is counted without 
multiplicity). If X is disjoint from its deck transformations 



C. 0{n) loop model 

Another model describing the statistics of loops on a lattice 
is the so-called 0(«) loop model^. Let us begin by defining 
the following 0(n) spin model on some (finite) lattice by the 
following partition function: 

r d^- 

Zoi„){x)= UT^W+^^rSj) (17) 

with S, G M", |S,| = 1 and £2„ is the «-dimensional solid an- 
gle. In order to obtain the loop model, we write S, • Sy = 

^I'^^y'^ + . . . + sl"^S^"^ and define n different colors (each of 
which will be associated with a specific component of the 
(9(«)-spins). Multiplying out all terms in Eq. ( fTTb . we have 
« choices for each bond plus a possibility of a vacant bond. 
Thus, the various terms are represented by an «-colored bond 
configuration; g — Qi,. . .Qn with denoting those bonds 

where the term sf'^S^P has been selected. Clearly, the various 
gi's are pairwise (bond) disjoint. Thus, for each g we obtain 
the weight 

W,;=Tr n xSl'^sfK.. n xS\"^S^;'\ (18) 

On the basis of elementary symmetry considerations it is clear 
that W(; 7^ if and only if each vertex houses an even number 
(which could be 0) of bonds of each color Once this con- 
straint is satisfied, we get an overall factor of x*'^' - with 
b{g) being the total number of bonds - times the product 
of the vertex factors obtained by performing the appropriate 
0{n) integrals. Obviously, these vertex factors depend only 
on how many different colors and how many of each of these 
colors enter each vertex (i.e. not on the particular colors in- 
volved nor on the directions of approach to the vertex). A 
particularly easy case is that of the honeycomb lattice where, 
due to a low coordination number, a maximum of two bonds 
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of a single color can visit a vertex. The corresponding vertex 
factor is then given by 



(19) 



leading (after summing over all n colors) to the following ex- 
pression for the partition function: 



{«} 



hia) 



(20) 



where ^(a) is the total number of occupied bonds (the total 
perimeter of all loops) while £(a) is the total number of loops. 
The last factor appears because each loop could be of one of 
n colors. The expression in ( |20l ) is well-defined for arbitrary 
n and x, so it can be taken as the definition of the 0{n) loop 
modeP^. 

We note that the partition function ( |20] | is exactly the norm 
of the ground state (^o |^o) of the quantum Hamiltonian // 
given by Eq. Q provided that x~n and (f- = n. 



D. Correlations 

Both of these models have a Coulomb gas representatior^ 
so that their correlation functions can be obtained from expo- 
nential operators in a Gaussian field theory with a background 
charge. Consider, for instance, the 0(«) model. Precisely the 
same loop expansion derived in section HV Cl ean be obtained 
from an SOS model on the kagome lattice. After integrat- 
ing out the triangular faces, the resulting SOS model for the 
heights on the hexagonal faces is has a low-temperature ex- 
pansion which is a sum over domain wall configurations. The 
weights of the SOS model (or, equivalently, 6-vertex model) 
are such that when a domain wall turns left, it acquires a fac- 
tor {x/n)e''^; when it turns right, a factor (x/n)e^'''. Since 
a loop will only close if the difference between the num- 
ber of right terms and the number of left turns is ±6, ev- 
ery closed loop receives a factor (x j n)^ e^^^~*- where h is the 
length of the loop. Summing over both orientations of the 
loop, we obtain 2(x/«)''cos6X- Hence, this is equivalent to 
the 0(n) loop model for /i < 2 if we take « = 2cos6X- When 
x> Xc = nl \/2 + \fi — «, this model is in its low-temperature 
phase, which is critical. 

The SOS model is a Coulomb gas with coupling g=-\ — 
6%/% and background charge —2(1 — g) (which ensures the 
correct phase factor for each turn of a loop). Consider the 
(Si -Sj) correlation function in the 0(n) model. Every con- 
figuration which gives a non-zero contribution must have one 
curve which does not close into a loop but has endpoints at 
spins at / and j. In Coulomb gas language, they have mag- 
netic charges ±1/2. In addition, they must each also have 
electric charge \ —g. Together with their magnetic charge, this 
will ensure that a factor e^^'"*- arises whenever the curve winds 
around either point (as required by the SOS vertex rules). This 
electric charge also cancels the background charge. Thus, the 
Coulomb gas operators corresponding to 0(«) spin operators 



have electric and magnetic charges (1 — g,±l/2). The ex- 
ponent associated with a correlation function between field 
with electric and magnetic charges (ei,mi) and (e2,»J2) is 
-e\e2/2g ~ gmim2/2. Hence, the 0(«) spin- 



^e^ ,mi ;e2 Jii2 



spin correlation function has the power-law decay^. 



(s('-)-s(o))^— , 



(21) 



where JCM = ig-^il-g)^- 

The loops which surround random clusters as described in 
Sec. lIV Al in the self-dual critical g'-state Potts model can also 
be mapped onto a 6-vertex model (and, therefore, a Coulomb 
gas) in a similar fashion. Precisely the same exponents are 
obtained. 

Although, they are crucial to our proof of gaplessness, these 
correlation functions are not the ones of direct physical in- 
terest. Physical correlation functions have a rather different 
behavior Consider equal-time correlations between quantum 
spins such as {o'-G~j) in the original quantum models ( I5I7I ). As 
we now argue, they are short-ranged in space. In the first 
model, the equal-time (c7; o~j) correlation function is related to 
the probability that a loop passes through ; and a loop which 
may or may not be distinct passes through j. Such correlation 
functions vanish in the 0{n) loop models. Analogously, in the 
second case, this correlation function is related to the proba- 
bility that the two corresponding bonds are parts of clusters, 
but not necessarily in the same cluster in the related q'-state 
Potts model. Such a correlation function once again vanishes. 

In Coulomb gas language for the associated statistical me- 
chanical models, the reason that such correlation functions 
vanish is that they are correlation functions of electrically neu- 
tral operators, such as gradients of the height (to which the 
local loop density corresponds). Such correlation functions 
vanish since they do not cancel the background charge. (The 
only exception is a height model, with central charge c = I, 
for which there is no background charge. Gradients of the 
height have power-law correlations.) As we have seen, alge- 
braic decay is possible for correlation functions of operators 
which are charged in the Coulomb gas picture, but these are 
non-local in terms of the spins o? since they measure, for in- 
stance, the probability that two spins o? and O'j lie on the same 

loop. (At d = 1,V2, this can also be seen from the fact that 
the ground state on the sphere is the same - and, therefore, 
has the same equal-time correlation functions - as that of a 
gapped Hamiltonia n"^-^" which is a sum of local commuting 
operators and, therefore, has correlation length zero.) 

Thus, the ground state wavefunction of Q has an under- 
lying power-law long-ranged structure which is apparent in 
its loop representation, but it is not manifested in the correla- 
tion functions of local operators a'^. As we will see momen- 
tarily, this long-range structure leads to gapless excitations 
for the Hamiltonian (|5]l and, therefore, long-ranged correla- 
tions in time in spite of the lack of long-ranged correlations 
in space. We call such a state of matter a quasi-topological 
critical point. 
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V. LOW-ENERGY EXCITATIONS 

In spite of the short-ranged nature of equal-time spin-spin 
correlation functions and the absence of any conservation laws 
for either of the Hamiltonians ( I5I7I ). we can construct a varia- 
tional argument that this general type of Hamiltonians is gap- 
less using the criticality of non-local correlation functions, 
specifically the scale-invariant nature of loops. 

The general idea of our proof is to produce a "twisted" state 
which is both orthogonal to the ground state and has a van- 
ishingly small expectation value of the Hamiltonian, in the 
general spirit of the Lieb-Schultz-Mattis theorem for quantum 
antiferromagnets^iiij^. 

For the purposes of this theorem we do not need to choose 
a particular Hamiltonian, but only require that it satisfy the 
following necessary conditions: 

• This is a c/-isotopy Hamiltonian whose ground state is 
described by Eq. ^ with a scale-invariant distribution 
of loops in the thermodynamic limit. We make the 
mathematically nontrivial assumption, widely accepted 
in physics, that the scaling limit exists. Furthermore, 
certain probability functions defined from the limit will 
be formally differentiated. These derivatives can easily 
be replaced by finite difference quotients, so continuity 
is, in fact, an adequate assumption. (We shall make this 
condition more precise later.) 

• A second important feature of the ground state which 
we use is that when viewed as a statistical mechanical 
ensemble by Prob(\|/, ) = it is, in the scaling limit, 
a critical system with the "no large cluster property," 
i.e. on an L X L square with periodic boundary con- 
ditions ("torus") VA, > 0,3e > so that the event "all 
loops within the random multi-loop \(// have breadth (as 
defined is Section ITVB I ) < XL" occurs with probability 
greater than e. (We have seen that the FK models for 
q> I have this property. We also expect this property 
to hold for the 0(n) loop models, 1 < « < 2; proving 
this would be quite interesting.) 

• The Hamiltonian is local: all allowed (i-isotopy moves 
have finite range (generalizing our theorem to the case 
of a quasi-local Hamiltonian whose terms decay ex- 
ponentially with increasing range is straightforward). 
Without loss of generality, we might as well assume that 
the range of all terms is limited to a single lattice pla- 
quette. 

• The terms responsible for the loop dynamics are 
bounded uniformly in the size of the system, i.e. for 
any two multiloops a and P, |(a|//f/|p)| < V. (We as- 
sume that the basis vectors |a) of the Hilbert space of 
Hd are orthonormal, i.e. (a|P) = unless a and P are 
identical multi-loops and (a|a) = 1). Notice that we 
make no assumptions about the other terms that might 
be there to enforce the multiloop constraint - these can 
potentially include hard-core interactions. 



There is little doubt that both of the presented Hamiltonians 
satisfy the above conditions. However, our argument is more 
complete in the second case (FK^, I <q <4). 

We will now construct our low energy excitation \\\ti) above 
the ground state |x|/o). We normalize so that (\(/o|i|/o) = 1, 
(i|/o|//|v|/o) = 0, and = 1. We will construct and 

then a family of "harmonics" all of norm one, and esti- 
mate the energy expectation values {\\tk\H\\\tk), A: > 1. In con- 
structing I'^k), we will use the language of the scaling limit 
for conceptual simplicity. It is quite routine, and we leave this 
to the reader, to back away from the scaling limit and write 
discrete formulae, replacing derivatives with difference quo- 
tients. 

We use C to denote a configuration (i.e. a multi-loop) near 
the scaling limit, i.e. L on T, the Lx L torus. By 

proposition |2] we know there is a function on configura- 
tions which remains almost surely defined in the scaling limit: 

= Our assumption is that, in the scaling limit, 

the probability p that a configuration C satisfied r^- < r is a 
continuous function p{r). The probability is computed with 
respect to the FK^j measure, i.e. with respect to |\|/o). As ex- 
plained earlier, we will formally treat p as differentiable, writ- 
ing dp/dr, but this may be treated as a difference quotient. 

The ground state wave function (on the trivial sector) |i|/o) 
is simply: 

|Vo)=Z-'/2£/(«)|a) (22) 
a 

where the sum is over all multi-loops a in the zero-winding- 
number sector and the normalization Z = L{a}'^^''"' is the 
partition function of the associated statistical mechanical 
model (Potts or 0(n)). 

We write the variational ansatz 

\Mfk) - Z-i/2£e27:«p('-a) /(a) (23) 
a 

where k is an integer. These states are orthonormal because 
a JO 

(24) 

In particular, \\\tit) with ky^Ois orthogonal to the ground state. 
Hence, the expectation value of the Hamiltonian in this state 
is an upper bound on the energy gap between the ground state 
and the first excited state. 

In the second line, we have used the fact that (a|i/rf|P) is non- 
zero only if £{a) = ^(P). Exchanging a with p exchanges each 
term in the sum in (|25] l with its complex conjugate. Hence, we 
can write: 

A < 2Re J i (e2'^'^'(P('-«)-''('"p)> - l) t/2^'(")(a|//,/|p) I 

I a,p;;p>roi J 
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Note that we have dropped the case rp = ra since this expres- 
sion vanishes for these configurations. The c/-isotopy Hamil- 
tonian can change the length of a loop by at most one lattice 
spacing, so (a|//^/|P) 7^ only if = ra + f where a is the 
lattice spacing. (Recall that ra has been defined in units of L.) 
Then, writing p{ra) - p{r^) « p'{ra) we have: 




Inserting /q drb{r — ra) = 1 into this expression, we have: 




where 

p(r) ^ -i ^ 8(r-ra)«f2«(«)(a|//,|P) 

a,P;rp>ra 

^ ^E8(r-r„)^2£(a)p_^ (27) 
^ a 

In the first line of dZTl ). we define p(r), which is the expecta- 
tion value of Pa, defined in the second line. If we can argue 
that p (r) < A for some constant A independent of r, then 

= 2aJ dr {\- cos, {2mkp'{r)a/L)) 

< 2A {I- cos {2KikMa/L)) (28) 

Here, we have taken < M, which follows from our mild 
continuity assumptions. Consequently, in the large L limit, 
A ~ /Lr. Note that the L^^ reproduces the classical scaling 
for the lowest eigenvalue of a string of length L, while the k^ 
dependence signals a quadratic dispersion relation associated 
to a soft mode. 

We are nearly finished. All that remains is to argue that p (r) 
defined by jZTl ) satisfies p(r) < A for some constant A inde- 
pendent of r. Pa is the number of distinct sites (i.e. distinct 
terms in H^i) where can produce a fluctuation stretching out 
the (x-direction of) the widest loop in the multi-loop a from 
breadth ra to breadth ^a + ^, i.e. by one lattice step. p(r) is 
the expectation value of pa, averaged over all configurations 
a whose widest loop has breadth r. Note that in the scaling 
limit p(r) can have no r dependence. Prior to reaching the 
scaling limit, we use argument r in p (r) to indicate maximum 
cluster breadth in units of -i. We now define pa.„ to be 1 if the 
widest loop K G a meets the right side of the unique smallest 
rectilinear box B, containing K, K C B, in n distinct "fingers" 
touching the right-hand wall and zero otherwise. We define 
p^(r) to be its expectation value, 

p;ir)^W5{r-ra)d^'^''^Pa,n (29) 
a 



Then 

00 

p(r) = Y.np;;{r) (30) 

n=l 

The Hamiltonian Hii can produce n distinct states contribut- 
ing to pj (r + |) for each state contributing to p„(r). For ex- 
ample, in the figure above, one may fluctuate any right-most 
finger further to the right. Hence, 

Pi('-+!)> £«P«W (31) 

n=l 

But the right-hand-side is simply p(r). Since Pj (r) is the en- 
semble average of p^^ j = 0, 1, it must satisfy < Pi{r) < 1. 
Hence, if there is a well-defined scaling limit, Pi(r) Pj, 
p(r) and p(r+ should converge. Then the latter is the de- 
sired constant A in the upper bound (l28b on the excited state 
energy since we can re-write dSTI ) as 

P < Pi (32) 

Naturally, this is only an upper bound, and in fact we have 
not proved these are not just other degenerate ground states. 
However, it seems unlikely that we have discovered ground 
state degeneracy (on the sphere) since there is no continuous 
symmetry of the Hamiltonian which could be broken or any 
other non-ergodicity. Furthermore, we expect co 0= to be 
the correct behavior of the low-lying excitations; our intuition 
is based on observations made for other quasi-topological 
models^i^MI, closely related to our d = 1 case. Due to this 
characteristic quadratic spectrum, such lines of critical points 
have been dubbed "quantum Lifshitz points'"*'^. 

Let us pause here and contemplate the physical reasons for 
having such gapless modes. These excitations are not Gold- 
stone bosons, since neither of the Hamiltonians (I5l7b pos- 
sesses any continuous symmetry (which could be broken). 
Rather, these gapless modes appear as a result of "bottle-neck" 
quantum dynamics which only allows loops to fluctuate by ei- 
ther slowly growing or slowly shrinking, one lattice plaquette 
at a time. As a result we can think of configuration space as a 
very elongated, essentially one-dimensional object parameter- 
ized by the diameter of the biggest loop. While the quantum 
dynamics is in principle ergodic, in order to reach a state with 
a long, order L loop from a state with only short loops, the 
entire "length" of this "worm-like" graph has to be traversed. 
From this analogy, we see that the eigenvalue problem for our 
Hamiltonian is very similar to the eigenvalue problem for the 
Laplacian operator on a string. (In contrast, geometries such 
as a complete graph or hypercube do have spectral gaps; their 
bonds tie them together more efficiently than links in a linear 
chain.) 

Now, what about the spectrum when the Jones- Wenzl pro- 
jectors are implemented^^? With the help of such additional 
terms, going from short to long loops can be achieved by 
merging existing loops together which can be done substan- 
tially faster than by "growing" them. These projectors would 
directly connect various points of the "worm-like" configu- 
ration space of the system. In the case d = I adding the two- 
strand projector essentially reproduces Kitaev's "toric code"— 
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and a gap is opened. (We expect this to be the case for c/ = — 1 
as well.) In the case of n — — 2, adding the three-strand 
projector is not sufficient because the probability of three long 
loops coming together within several lattice constants from 
each other - a necessary condition for a JW projector to act 
efficiently - vanishes as L^" for some a > 0^. Hence the ef- 
fect of such terms on the spectrum is too weak to open a gap 
(although it might lead to the "stiffening" of gapless excita- 
tions by reducing the dynamical exponent z). 

VI. CONCLUSIONS 

In this paper, we have analyzed the spectrum of low lying 
excitations for a general class of local Hamiltonians whose 
ground state(s) are characterized by t/-isotopy. Using the sta- 
tistical properties of their non-local degrees of freedom we 
established the analog of a Lieb-Schultz-Mattis theorem for 
quasi-topological systems. The excitations are gapless, with 
the variational ansatz strongly suggesting a quadratic (CO °^ k^) 
dispersion. What may strike one as interesting and counter- 
intuitive is the fact that both the Hamiltonian and the corre- 
lations of local operators are perfectly short-ranged, however, 
the quantum dynamics is constrained to operate on non-local, 
quasi-long-ranged objects - loops - which in turn leads to a 



gapless spectrum. Finally, it is also interesting to remark on 
the potential implications of such a behavior from the per- 
spective of understanding quantum glasses. The idea of using 
similar quasi-topological models to model glassy behavior is 
not entirely new'*'' but until very recently^, all proposed 
models had a serious drawback, namely quasi-long range cor- 
relations between local degrees of freedom. This is contrary 
to the fact that experimentally observed slow glassy dynamics 
has not been accompanied by any divergent correlations. In 
this paper we have explicitly demonstrated that such behavior 
is entirely possible in the context of quasi-topological quan- 
tum critical points. 
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A multiloop is a collection of non-intersecting loops and arcs on a 
surface. The end points of an arc are required to lie on the bound- 
ary of the surface, see Ref.[l2l. 

See Ref. |4^ where private communication with A. Kitaev and 
(psKuperberg is referenced. The required family of commuting 
projectors is easily derived in the Turaev-Viro approaches by writ- 
ing surface X interval, Z x / = handle-body union 2-handles. The 
disjoint attaching curves of the 2-handles yield the commuting 
projectors. 

Here we use the term "frustration" in reference to a Hamiltonian 
which can be written as a sum of projectors yet has no zero- 
modes. In this sense, new terms breaking the extensive degen- 
eracy down to a crystalline structure will introduce frustration. 
We stress that this is a term-by-term correspondence in the sense 
that, given the appropriate basis, the squared amplitudes of the 
basis states comprising the ground state are identical to the Gibbs 
weights of the corresponding statistical mechanical states (up to 
an overall constant). 

In a fully-packed case, there is an additional geometric constraint 
limiting these winding numbers to either only even or only odd, 
depending on the size of the torus 

This derivation can be made slightly more formal if one defines a 
joint measure on both spin and bond configurations: p{a,(£)) = 
v''(")a(<t, to). This is known as Edwards-Sokal measured. As 
follows from Eq. l ll2b . when traced over both spins and bond oc- 
cupations, it gives the desired partition function. Its marginal with 
respect to bond configurations is the Gibbs weight for spins, while 
its marginal with respect to spin configurations defines the weight 
for random clusters. 

It is also known as a cyclomatic number and for the case of a 
planar graph is equal to the number of finite face^. 



